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We study coherent quantum transport through a superconducting film connected to normal-metal 
electrodes. Simple expressions for the differential conductance and the local density of states are 
obtained in the clean limit and for transparent interfaces. Quasiparticle interference causes periodic 
vanishing of the Andreev reflection at the energies of geometrical resonances, subgap transport, and 
gapless superconductivity near the interfaces. Application of the results to spectroscopic measure- 
ments of the superconducting gap and the Fermi velocity is analyzed. 



I. INTRODUCTION 



Conducting properties of normal metal/superconductor/normal metal (NSN) heterostructures have attracted con- 
siderable interest^*^ Electronic transport through a normal metal/superconductor (NS) junction, with an insulating 
barrier of arbitrary strength at the interface, was studied by Blonder, Tinkham, and Klapwijk (BTK),'^ and the 
Andreev reflection^ was recognized as the mechanism of normal-to-supercurrent conversion^i^ The BTK model can 
be used when electrons pass incoherently through the interfaces of NSN double junctionsii 

However, in clean superconducting heterostructures coherent quantum transport is strongly influenced by size 
effects due to the interplay between geometrical resonances and the Andreev reflectioni2i2iiii Well-known examples 
are the current-carrying Andreev bound statesii*i^ and multiple Andreev reflectionsi^iMii^ in superconductor/normal 
metal/superconductor (SNS) junctions. Since early experiments by Tomaschji& the geometric resonance nature of the 
differential conductance oscillations in SNS and NSN tunnel junctions has been attributed to the electron interference 
in the central layer jiiiSiiSi^'^ Recently, McMillan- Rowell oscillations were observed in SNS tunnel junctions of d-wave 
superconductors, and used in measurements of the superconducting gap and the Fermi velocity^ 

In this paper, we focus on size and coherency effects in NSN double junctions with a clean conventional (s-wave) 
superconductor and transparent interfaces. In that case, instead of a fully numerical treatment jSiSS a simple analytic 
form for the differential conductance and the local quasiparticle density of states (DOS) can be obtained from the 
solution of the Bogoliubov-de Gennes equation. The proximity effect, which is significant in thin superconducting 
films, is taken into account through the self-consistent calculation of the pair potential. Striking consequence of the 
quasiparticle interference in the clean metallic NSN junctions include periodic vanishing of the Andreev reflection at 
the energies of geometrical resonances and the characteristic changes of DOS with respect to the bulk superconductor. 



II. SCATTERING PROBLEM 



We consider an NSN double junction consisting of a clean superconducting layer of thickness d, connected to 
normal-metal electrodes by transparent interfaces. The quasiparticle propagation is described by the Bogoliubov-de 
Gennes equation 

Fo(r) A(r) \ /^r)\ _ r,(u{v)\ 
A*(r) -i7o(r) ) \v{v) ) " \v{v) 

with Ho{r) = —K^V^ /2m — Ep, where Ep is the Fermi energy. For simplicity, the magnitude of the Fermi wave 
vector, kp = \/2mEp /h? , is assumed to be constant through the junction. The superconducting pair potential is 
taken in the form A(r) = A9(z)9(Z — z), where 0(z) is the Heaviside step function, and the z-axis is perpendicular 
to the layers. In Eq. E \s the quasiparticle energy with respect to the Fermi level. The electron effective mass m 
is assumed to be the same for the entire junction. The parallel component of the wave vector k|| is conserved, and 
the wave function 



u{y) 
v{v) 



exp(ik[| • r) V(^), (2) 



satisfles the continuity boundary conditions for ili{z) and tp'{z) at z = and z ~ d. The four independent solutions 
of Eq. ^ correspond to the four types of injection: an electron or a hole, from cither the left or the right electrodei^ 
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For the injection of an electron from the left, with energy E > and angle of incidence 9 (measured from the 
z-axis), solution for ip{z) has the following form 

exp(zA;+z) ^ q ^ + 0) exp{ik~z) ^ ]' ^ z < 0, 

iPiz) ^ { bi{E,0)exp{iq+z) (^^^ +b2{E,0)exp{iq'z) (^^^ < z < d, (3) 

c(i?, 0) exp(ifc+z) ^ Q ^ z > d. 

Here, u = ^^(1 + n/E)/2 and v = ^(1 - n/E)/2 are the BCS coherence factors, and n ^ The z- 

components of the wave vectors are = ^ {2m/h?){Ep ± E) — k^| and = ^ {2m/h?){Ep ± fi) — , where 

|k|| I = ^y {2m/h?){Ep + E) smO. The coefficients a and c are, respectively, the probability amplitudes of the Andreev 
reflection as a hole of the opposite spin (AR) and transmission to the right electrode as an electron (TE). The 
amplitudes of the Bogoliubov electron-like and hole- like quasiparticles, propagating in the superconducting layer, are 
given by the coefficients hi and 62- The normal reflection of electrons and the transmission to the right electrode as a 
hole are absent. Solutions for the general case of insulating interfaces and the Fermi velocity mismatch are given in 
Ref. 7. 

By neglecting the small terms E/Ep <C 1 and A/Ep <C 1 in the wave vectors, except in the exponent 

C = d{q+-q-), (4) 

solutions of the boundary-condition equations can be written in a simple form: a{E,9) = uv[l — exp{iCz/d)]/[u^ — 
v'^exp{iCz/d)], bi{E,9) = u/[u'^-v^ exp{iCz/d)], b2{E,9) = -v exp{iCz/d)/[u'^ -v^ exp{iCz/d)], cindc{E,9) = -i{v?- 
v'^) ey.Y>{iC,z / d) /[v?' — exp(iCz/c?)]. From the probability current conservation, the probabilities of outgoing particles 
satisfy the normalization condition A{E,e) + C{E,9) ^ 1, where A{E,9) = \a{E,9)\'^ and C{E,9) = \c{E,9)\'^. 
Explicitly, the AR probability is 



AiE,9) 



Asin(C/2) 



£:sin(C/2) -h jl7cos(C/2) 



(5) 



Solutions for the other three types of injection can be obtained by the same procedure. In particular, if a hole 
with energy —E and angle of incidence 9 is injected from the left, the substitution g"*" ^ q~ holds, and the scattering 
probabilities are the same as for the injection of an electron with E and 9. Therefore, in order to include the description 
of both electron and hole injection, the calculated probabilities should be regarded as even functions of E. Also, for 
an electron or a hole injected from the right, the probabilities are the same as for the injection from the left. 

From Eq. ((SJ, it follows that A{E, 9) = 0, and consequently C{E, 9) = 1, when 

C = 2n7r (6) 

for n = 0, ±1,±2, . . .. Therefore, the Andreev reflection vanishes, and the electron transmission becomes complete 
(without creation or annihilation of Cooper pairs), at the energies of geometrical resonances in the quasiparticle 
spectrum, given by Eq. ©. 

For thin superconducting films, the proximity effect is important and the self-consistent treatment of the pair 
potential is required. However, the spatial variation of the pair potential is negligible for c?/^o ^ 1, whereto = hvp/nAo 
is the BCS coherence length in the bulk superconductor, and we can use the previous solution with the step-like form 
for A. The pair potential is given by the self-consistency equation 

d^'k\\ j dVL u{Y)v*{v)iimk{El2kBT), (7) 

where A is the BCS coupling constant and iV(0) = mkp /2Tr'^h'^ is the normal-metal density of states (per spin) at the 
Fermi level. We calculate the spatial average of A{z) using the standard iteration procedure 

1 



A^+l - - / A,{z)dz, (8) 



by setting fl = ^ E'^ — A? in Eq. Q in order to obtain Ai{z) in the i-th iteration. ^'^ Starting from the bulk value 
Aq, we repeat this procedure until the difference between A^+i and A^ becomes sufficiently small. In Fig. ^ we show 
the resulting A as a function of d, at zero temperature. 
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III. DIFFERENTIAL CONDUCTANCE AND DENSITY OF STATES 

The normalized differential conductance at zero temperature for a planar double junction is 



/ d0 sin cos [1 + A{E, 0)] , 
Jo 



Gn 

where Gn = e^//i. In the one-dimensional (ID) case, the normalized differential conductance is simply 

9iE) 



(9) 



G 



N 



= 1 + A{E,0). (10) 



Coherent transport through an NSN double junction is influenced by considerable changes of the quasiparticle 
density of states with respect to the bulk superconductoriiii^ Adopting the method introduced by McMillan,— 
Ishii^ and Furusaki and Tsukada,^ we evaluate the Green functions by combining the solutions of the Bogoliubov-de 
Gennes equation. From the imaginary part of the retarded Green function, we obtain the local value of the partial 
density of states (PDOS) for the superconducting film, normalized with respect to the ID normal-metal density of 
states per spin at the Fermi level. 



N{z,0,E)^-—^ Re \2EHE^ + n^) + 2E^A''cosC 

i (E) COS0 I 

+ [cos{C{z/d - 1)) + cos(Cz/d)][A'^ - A^{E^ + n"^) cosC] 
+2E^A^[sm{C{z/d~l))~sm{Cz/d)]smC^ , (11) 

where 

T{E) = [{E^ + r^^) cos C - A^]^ + AE^n^ sin^ (. (12) 
For a planar junction, the local DOS for the superconducting film is given by 

N{z,E) _ 



7V(0) 



/ d0 sm0 COS0 N{z,0,E). (13) 
Jo 



Here, N{z, 0, E) — N{z, E)/N{0) — 1 in the normal-metal electrodes {z < and z > d). 

Characteristic features of the coherence in NSN double junctions are subgap transport of electrons and oscillations 
of the conductance. Due to the interference effect, the Andreev reflection is suppressed for E < A, whereas the AR 
probability oscillates with E and d, for E > A. The subgap transport is dominant in thin superconducting fllms, 
and the oscillatory behavior is apparent in thick films. These oscillations are more pronounced in the ID case, which 
is shown in Fig. |21 It can be clearly seen that the positions of the minima in g{E) exactly match the positions of 
the maxima in PDOS, corresponding to the geometrical resonances imposed by Eq. J^. Results can be applied to 
spectroscopic measurements of A and vp. From Eq. (jSJ for 9 = 0, the energies of the geometrical resonances. En, 
satisfy 



where n = 1,2, . . . counts the conductance minima. Therefore, the plot of E^ vs has the intercept equal to A^, 
and the slope equal to tan~^ [(Trftwi^/d)^] . An example is shown in the Inset of Fig. |21 Note that the points obtained 
from the minima of G{E) follow practically the same linear law as those of g{E). 

The obtained results are in a strong contrast with the conductance spectrum and DOS of tunnel junctions. In the 
latter case, g{E) and PDOS reduce to (5-function-like spikes at the bound state energies of the isolated superconducting 
film, given by Iq'^ = uitt and lq~ = ri27r, where ni — n2 = 2n, with n given in Eq. ©.^'^ 

The conductance spectra G{E) /Gn and the local DOS for the superconducting film in a planar NSN double junction 
are illustrated in Figs. ElandQl for d/^o = 1 and d/^o = 10, using the self-consistent pair potential A/Ao = 0.816 and 
A/Aq = 1, respectively. In Fig.O the conductance spectra are compared to the BTK result describing the incoherent 
transport. It can be seen in Fig. ^a) that the energy spectrum of a thin superconducting film is practically gapless. 
For a thick superconducting film, this is the case only near the interfaces. Fig. 0Jb). The characteristic oscillations 
of the local DOS become more apparent for a thick film. The oscillations of G{E) and of the local DOS are not 
completely smeared in the 3D case. The positions of the minima in G{E) or the maxima of the local DOS are still 
close to the positions of the resonant quasiparticle states in the ID case. 
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IV. CONCLUSION 



We have derived simple expressions describing the finite size and coherency effects in the metallic NSN double 
junctions, suitable for experimental data analysis. The main features of the coherent quantum transport through a 
superconducting film are subgap transport and oscillations of the differential conductance. The Andreev reflection is 
suppressed below the gap, especially in thin films. The conductance oscillates as a function of the layer thickness and of 
the quasiparticle energy above the superconducting gap. Periodic vanishing of the Andreev reflection at the energies 
of geometrical resonances is an important consequence of the quasiparticle interference; the conductance minima 
correspond to the maxima in the density of states. Measurements of the conductance spectra of metallic NSN double 
junctions could be used for accurate determination of the pair potential and the Fermi velocity in superconductors. 

We are grateful to Ivan Bozovic for useful discussions. This work has been supported by the Serbian Ministry of 
Science, Technology, and Development, Project N°1899. 
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FIG. 1: The pair potential A as a function of the film thickness d, at zero temperature. 

FIG. 2: (a) Partial {6 = 0) conductance spectra, g{E)/GN, and (b) corresponding PDOS, N{d/2,0,E), for d/^o = 10. Inset: 
square of resonant energies, E^, vs obtained from the minima of g{E) (open circles) and G{E) (full squares). 

FIG. 3: Differential conductance spectra G{E)/Gn for d/^o = 1 and 10 (solid curves). The BTK result is shown for comparison 
(dashed curve). 

FIG. 4: Local DOS, N{z,E)/N{0), for z — d/2 (solid curves) and z = or d (dotted curves) as a function of E/Ao is shown 
for d/Co = 1 (a), and d/^o = 10 (b). 



